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(57) Abstract: The invention concerns a device and a method 
assisting in^monitoring andVor simulating a complex system, in 
particular a telecommunication network. It consists in storing 
data in matrix form, of dynamically variable structure, with it- 
erative multiplication in max -plus algebra of a current matrix of 
network data and an instantaneous matrix of the network param- 
eters. 

(57) Abreg£ : L' invention concerne un dispositif et un proce'de' 
d'aide a la surveillance et/ou de simulation d'un systeme com- 
plexe, notamment d'un reseau de communication. On prevoit un 
stockage des donnees sous forme matricielle, de structure dyna- 
miquement variable, avec une multiplication iterative dans l'al- 
gebre MAX-PLUS d'une matrice courante de donnees de reseau 
et une matrice instantanee de parametres du reseau. 
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surveillance et simulation perf ectionn6es de systfemes 

nnm plPxes, notamment de m^canismes e t de controles de flux et 
5 de congestions dans des r^seaux de communication 

L' invention concerne la surveillance et la simulation de 
systfemes complexes . 

10 Dans le cadre du controle de flux et de congestions dans des 

r^seaux de communication , notamment de type internet, une 

analyse fine du d<§bit offert est souhait<§e pour estimer les 
influences respectives des parametres du r6seau* 

15 Avec le d£veloppement des techniques de communications en 
ligne et les problemes rencontres, notamment de congestion, 
diff^rents protocoles de contrdle de flux et de congestion 
ont vu le jour, notamment le controle TCP (de 1* anglais 
"transmission control protocole" ) • 

20 

On connait des proc6d6s d' analyse de ces protocoles- Parmi 
ces proc6d6s connus, un proc£d£ bas6 sur une expression 
math«§matique adapt<§e du d6bit dans un protocole de type TCP 
a permis une approche analytique du controle. Le principe sur 
25 lequel repose ce procede est decrit notamment dans z 

Mo Mathis, J- Semske, J. Mahdavi and T. Ott, "The 
Macroscopic Behavior of the TCP Congestion Avoidance 
Algorithm", Computer Communication Review, 27(3), July 
30 (1997)o 

Ce proc6d6, quoique prometteur, a montr6 ses limites dans des 
applications pratiques, notamment le fait que le caract&re 
al^atoire du trafic ne soit que partiellement pris en compte, 
35 ou encore qu'il n6cessite une approximation de tous les 
noeuds du r^seau h un unique noeud Equivalent, virtuel. 

Un autre procEdE connu, prenant mieux en compte 1* approche 
stochastique, a permis d ' apprEhender le caractfere al6atoire 
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du trafic. Ce proc<§de, plus recent, d^coule du principe 
d^crit dans les travaux suivants : 

Padhye, J., Firiou, V. , Towsley, D. and Kurose, J. 
5 "Modeling TCP throughput: a simple model and its empi- 

rical validation", Proc. of ACM SIGCOMM (1998) 

Padhye J., Firiou V. , Towsley D., "A Stochastic Model of 
TCP Reno Congestion Avoidance and Control", Technical 
1 q Report, 99-02, CMPSCI, Univ. of Massachusetts, Amherst 

(1999) . 

Cependant, ce proc<§de a aussi montr€ ses limites, notamment 
le fait qu'il necessite encore une approximation de tous les 
15 noeuds du reseau a un noeud Equivalent. 

La pr^sente invention vient am£liorer la situation. 

Selon une approche diff^rente, 1' invention propose d'utiliser 
20 une representation dans 1 9 algebre dite "max-plus" de systfemes 
complexes, tels que des r^seaux de communication et notamment 
du controle de flux et de congestions. 

Pour obtenir le detail des principes math^matiques sur 
25 lesquels repose une telle representation, on pourra se 
ref^rer a l'ouvrage suivant : 

F. Baccelli, G. Cohen, G.J. Olsder, and J. P. Quadrat, 
Synchronization and Linearity, Wiley (1992). 

30 

Globalement, 1' algebre max-plus scalaire est un semi-anneau 
sur la droite reelle ou 1' addition devient la fonction "max" 
(valeur la plus grande parmi un ensemble de valeur) et la 
multiplication, la fonction "plus" (some). L' utilisation de 
35 1' algebre max-plus permet de ramener les calculs d'un systeme 
compliqu6 a une simple representation matricielle. 

La Demanderesse a montr6 et verifi£ en pratique que 1' utili- 
sation de 1* algebre max-plus s ' adapte de fa?on tres satis- 
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faisante h la surveillance et a la simulation de systemes 
tels qu'un reseau de communications, control^ ou non. Elle 
permet notamment de surmonter le caractere aieatoire des 
parametres du r6seau # tout en considerant une plurality de 
noeudSo De plus, la Demanderesse a montre que la representa- 
tion d'un reseau utilisant un protocole TCP etait lineaire 
dans 1'algSbre max-plus , ce qui permet, en pratique, d'appli- 
quer des traitements de donnees simples. 

La presente invention porte alors sur un dispositif d'aide & 
la surveillance et/ou de simulation d'un systeme complexe, 
notamment d'un reseau de communication. 

Selon une premiere importante de 1' invention, le dispositif 
comprend s 

- une m^moire pour stocker de premieres donnees representa- 
tives de parametres du reseau, ainsi que pour recevoir au 
moins de secondes donnees representatives d ' ev^nements dans 
le reseau, une portion de ladite memoire etant r6serv6e au 
stockage de donnees sous forme matricielle, 

- un module de calcul, apte & effectuer sur au moins deux 
matrices de structure dynamiquement variable, une operation 
formant produit selori 1 ' algebre dite MAX— PLUS , 

- un module de modeiisation pour construire au moins une 
premiere matrice et une matrice courante respectivement en 
fonction des premieres donnees et des secondes donnees, selon 
un module choisi, et 

- un module pilote pour appliquer repetitivement la premiere 
matrice et la matrice courante au module de calcul, la 
matrice produit obtenue devenant une nouvelle premiere 
matrice <, 

Pr6f erentiellement , les premieres donnees comprennent des 
informations sur la topologie du reseau, telles que le nombre 
de routeurs traverses par la connexion ci surveiller ou & 
simuler, les proprietes de ces routeurs (tailles memoires 
dites "buffers", ou autre), les proprietes statistiques des 
trafics offerts dans le reseau, etc. 
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Selon une autre caracteristique avantageuse de 1' invention , 
le module de moduli sat ion comprend : 

- un sous-module de modelisation statique, pour construire la 
premiere matrice en fonction des premieres donnees, et 
5 - un sous-module de modelisation dynamique, pour construire 
au moins une matrice courante en fonction des secondes 
donnees . 

Avantageusement, le dispositif de 1 ' invention est apte & 
10 traiter des matrices comportant des coefficients dynamique- 
ment variables, dont au moins la matrice courante pr6citee. 
Dans le traitement qu'effectue le dispositif selon 1' inven- 
tion , les matrices construites sont avantageusement de m§me 
dimension. 

15 

Ces secondes donnees comprennent pref erentiellement des 
informations relatives k des pertes dans le reseau, a des 
flux transverses dans le reseau, par rapport a une connexion 
contr6iee ^ surveiller ou a simuler, a des congestions dans 
20 le reseau, ou encore a des depassements de deiais dans le 
reseau dits "time-out". 

Selon une autre caracteristique pref erentielle avantageuse , 
la matrice produit obtenue est un vecteur represente par une 
25 matrice a colonne unique, ce qui permet de limiter les 
traitements et leur duree. La premiere matrice, representa- 
tive des parametres du reseau, est avantageusement structure 
au depart comme un vecteur. 

30 Selon une seconde caracteristique importante de 1* invention, 
la matrice produit obtenue est representative d'un debit dans 
le reseau associe ci la connexion a surveiller ou & simuler, 
d'un d^bit moyen dans le reseau, ou encore de fluctuations 
d'un debit instantane dans le reseau. 

35 

Selon une troisieme caracteristique importante de 1 ' inven- 
tion, le module de modelisation est agence pour construire 
successivement une pluralite de matrices, en norabre corres- 
pondent sensiblement au nombre de paquets dans le reseau. 
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Le modele choisi comprend pr<§f §rentiellement la consideration 
de la taille variable d'une fenetre utilis^e pour controler 
le nombre de paquets dans le r6seau. 

5 II peut s'agir d'un r6seau control^ par un protocole de type 
TCP, comprenant typiquement des routeurs k discipline de type 
"premier arrive premier servi", ou encore des routeurs a 
discipline de type WFQ (de 1' anglais "weighted fair queu- 
ing") o Le protocole TCP controlant le r£seau peut aussi bien 
10 etre base sur un module de Reno ou un modele de Tahoe, comme 
on le verra plus loin, 

Le service du r£seau peut etre deterministe, ou encore 
aieatoire, comme on le verra en detail plus loin- 

15 

La pr^sente invention vise aussi un procede d'aide h la 
surveillance d'un systeme complexe, notamment d'un r6seau de 
communication. Un tel precede comprend globalement les Stapes 
suivantes s 

20 a) obtenir des premieres donnees representatives de paramfe- 
tres du reseau, 

b) construire une premiere matrice / selon un module choisi, 
en fonction desdites premieres donnees, 

c) recevoir, & un instant choisi, au moins de secondes 
25 donnees representatives d'evenements dans le reseau, 

d) construire au moins une seconde matrice de structure 
dynamiquement variable, selon le module choisi, en fonction 
des secondes donnees, et 

e) effectuer sur lesdites matrices une operation formant 
30 produit selon l'algfebre dite MAX-PLUS, la matrice produit 

obtenue etant representative de 1 ' 6tat du reseau audit 
instant choisi. 

S'il est souhaite de suivre une evolution temporelle de 
35 l'etat du reseau £ des instants choisis, le procede comporte 
avantageusement l'etape suppiementaire suivante s 

f) repeter, & des instants choisis, les Stapes c), d) et e), 
tandis que la matrice produit obtenue devient la premiere 
matrice aprfes l'etape e). 
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La presente invention vise aussi un proc6de de simulation 
d'un systeme complexe, notainment de nuScanismes et de contro- 
les de f lux et de congestions dans un rfiseau de communica- 
tion. Ce procede comprend globalement les Stapes suivantes : 
5 a) obtenir des premieres donn^es representatives de parame- 
tres propres au reseau, 

b) construire une premiere matrice, selon un modele choisi, 
en fonction desdites premieres donn^es, 

c) simuler des evenements dans le reseau et pr6voir au moins 
10 de secondes donn^es representatives desdits Evenements , 

d) construire au moins une seconde matrice selon le module 
choisi, en fonction desdites secondes donn^es, et 

e) effectuer sur lesdites matrices une operation formant 
produit selon 1'algebre dite MAX-PLUS, la matrice produit 

15 obtenue etant representative d'un 6tat du reseau subissant 
lesdits evenements. 

Pour prevoir une evolution de l'etat du reseau en fonction 
des evenements qu'il subit, ce procede comporte avantageuse- 
2 0 ment l'^tape supplementaire suivante : 

f) r<§p£ter, pour des evenements successifs, les Stapes c) , d) 
et e ) , tandis que la matrice produit obtenue devient la 
premiere matrice apres l'etape e) . 

25 D'autres caracteristiques et avantages de 1' invention 
apparaitront a 1'examen de la description detaillee ci-apres, 
et des dessins annexes sur lesquels : 

- la figure 1A repr^sente schematiquement un dispositif au 
30 sens de la presente invention, 

- la figure 1 repr^sente schematiquement un nombre K de files 
en tandem dans un reseau, avec controle de flux, 

35 - la figure 2A represente un exemple d' evolution pas-a-pas de 
dateurs et de la taille d'une fenetre utilisee pour contrdler 
le nombre de paquets dans le reseau, 



iNSDOCID: <WO 0165772A1 l_> 



WO 01/65772 PCT/FR01/00579 

7 

- la figure 2 represente des interactions entre plusieurs 
paquets dans le reseau, 

- la figure 3 represente une variation du debit (courbe en 
5 traits pleins ) , obtenue par simulation, dans un reseau de 

protocole TCP base sur le modele de Tahoe sans la phase 
exponentielle , 

- la figure 4 illustre une interpretation graphique des 
10 debits asymptotiques dans un reseau de protocole TCP base sur 

le module de Reno avec deiais deterministes , 

- la figure 5 represente une variation du debit, obtenue par 
simulation et montrant une decroissance du debit en cas de 

15 pertes aieatoires dans le reseau, 

- la figure 6 repr6sente une variation du debit (courbe en 
traits pleins), obtenue par simulation, dans un reseau de 
protocole TCP base sur un modele de Reno markovien, 

20 

- la figure 7 represente des variations de debit comparees, 
obtenue par simulation, dans des reseaux de protocoles TCP 
bas^s respectivement sur un modele de Reno deterministe RD 
(traits pleins), de Reno markovien RM (traits pointings 

25 longs), de Tahoe deterministe TD (traits pointilies moyens ) 
et de Tahoe markovien TM (traits pointings courts), 

- la figure 8 represente une variation du debit (courbe en 
traits pleins), obtenue par simulation, dans un reseau de 

30 protocole TCP base sur un module de Tahoe avec phase expo- 
nentielle, 

- la figure 9 represente des variations comparees de debits, 
obtenues par simulation, dans un reseau de protocole TCP base 

35 sur un module de Tahoe avec des services s 3 et s 8 respective 
ment constants et egaux a 1, et 
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- la figure 10 repr^sente schematiquement un r^seau avec ses 
files et ses routeurs . 

L ' annexe I comprend les formules et equations El E28 
5 auxquelles la description d^taillee ci-apres fait reference. 

L' annexe II comprend les references bibliographiques [1] & 
[13] indexees entre crochets dans la description ci-apres. 

10 Les dessins, les annexes et la description ci-apres contien 
nent, pour 1'essentiel, des elements de caractere certain, 
lis pourront done non seulement servir a mieux f aire compren- 
dre la pr^sente invention, mais aussi contribuer a sa 
definition, le cas gcheant . 

15 

En se r^ferant a la figure 1, le dispositif se presente sous 
la forme d'un ordinateur comprenant une unit6 centrale UC 
munie d'un microprocesseur pP qui coop^re avec une carte-mere 
CM, Cette carte-mere est relive a divers 6quipements, tels 

20 qu'une interface de communication COM (de type Modem ou 
autre), une memoire morte ROM et une memoire de travail RAM 
(memoire vive). La carte-mere CM est reliee en outre a une 
interface graphique IG, laquelle pilote l'affichage de 
donn6es sur un ecran ECR que comporte le dispositif. II est 

25 prevu en outre des moyens de saisie, tels qu'un clavier CLA 
et/ou un organe de saisie dit "souris" SOU, relics & l'unit^ 
centrale UC et permettant a un utilisateur une interactivity 
avec le dispositif. 

30 La memoire ROM, ou encore la memoire RAM stocke les premieres 
donn^es pr£cit6es, representatives des parametres du reseau 
(topologie, propri^tes des routeurs, etc). Dans l'exemple, la 
memoire RAM regoit les secondes donn^es pr£citees, represen- 
tatives d'6v6nements dans le reseau (flux transverses, 

35 congestions, pertes, etc). Dans le cadre d'une aide h la 
surveillance du reseau ces secondes donn^es peuvent etre 
regues par 1' interface de communication COM. Dans le cadre 
d'une simulation, 1 ' acquisition de ces secondes donn^es peut 
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etre effective par un calcul base sur un module de simula- 
tion, comme on le verra plus loin. 

La memoire RAM au moins peut etre adressable en fonction de 
5 rang<§es et de colonnes de matrices et permettre ainsi un 
stockage de donnees sous forme matricielle. 

La m^moire ROM comprend un module de moduli sat ion MOD qui, en 
cooperation avec le microprocesseur yP, permet de construire 
10 la premiere matrice pr^citee et une matrice courante, 
respectivement en fonction des premieres donnees et des 
secondes donnees, selon un modele choisi que l'on verra plus 
loin. 

15 La memoire ROM comprend un module CAL qui, en cooperation 
avec le microprocesseur pP, permet d'effectuer sur au moins 
deux matrices de structure dynamiquement variable, une 
operation formant produit selon 1 ' algebre MAX-PLUS - Dans 
l'exemple, on entend par "matrices de structure dynamiquement 

2 0 variable", des matrices dont les coefficients au moins sont 
dynamiquement variables- Avantageusement , les modules qui 
seront decrits ci-apres permettent de ramener les matrices 
construites (et plus particulierement les matrices courantes) 
& des matrices dont seuls les coefficients sont dynamiquement 

25 variables o 

Le module de modelisation MOD comprend alors s 
- un sous-module de modelisation statique ST, pour construire 
la premiere matrice en fonction des premieres donnees, et 
30 - un sous-module de modelisation dynamique DYN, pour cons- 
truire, en fonction des secondes donnees, au moins une 
matrice courante dont les coefficients sont dynamiquement 
variables o 

35 La memoire ROM comprend en outre un module PIL qui, en 
cooperation avec le microprocesseur pP, permet d'appliquer 
repetitivement la premiere matrice precitee (comprenant les 
premieres donnees) et la matrice courante (comprenant les 
secondes donnees) au module de calcul CAL. La matrice produit 
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obtenue est stockee «n m6moire et devient une nouvelle 
premiere matrice. Elle peut etre. ensuite multipliee (dans 
l'algebre MAX-PLUS) a une autre matrice courante, comportant 
de nouvelles secondes donnees representatives de nouveaux 
evenements dans le reseau. 

Diverses approches ont ete proposees pour apprehender les 
proprietes cles du mecanisme de controle de flux & fenetre du 
type TCP a partir notamment de considerations heuristiques , 
de simulations, d ' approximations fluides ou encore d' analyses 
markoviennes references [10,11,1,12,13/14]. Tous les modeles 
analytiques sont bases sur la reduction du reseau a un seul 
noeud representant le goulot d ' etr anglement [9]. 

Par ailleurs, il a ete recemment demontr6 que le controle de 
flux a fenetre d'un reseau multidimensionnel admet une 
representation max-plus lineaire lorsque la taille de la 
fenetre est constante (reference [5])- Ici, le Demandeur 
s'est interesse pref erentiellement aux modeles qui combinent 
le mecanisme de controle adaptatif de TCP et un reseau 
multidimensionnel constitue de plusieurs routeurs en serie. 
La dynamique d'un tel xeseau controle est avantageusement 
decrite au niveau "paquets" via des iterations de produits 
matriciels dans l'algebre max-plus. On considere ici a la 
fois le cas ou les temps de transmission des paquets sont 
deterministes et les divers modeles stochastiques qui ont ete 
utilises dans l'etat de la technique, en particulier, les cas 
ou il y a des pertes aieatoires en plus des pertes dues au 
depassement de la capacite des memoires tampons (ou "buf- 
fers"), et avantageusement le cas ou les temps de transmis- 
sion des paquets sont aleatoirement perturbes par les autres 
traf ics . 

Beaucoup d' aspects cles du protocole peuvent etre represen- 
tes: pertes de congestion, pertes aieatoires, deiais (ou 
"time-outs") de propagation ou deiais dus aux attentes ou 
encore au mecanisme de controle de flux, etc* Comme on le 
verra plus loin, cette approche permet d'obtenir des formules 
explicites pour le debit maximum alloue, lorsque les pertur- 
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bations sont deterministes ou aleatoires . Ces f ormules sont 
asymptotiquement compatibles avec les formulas connues 
lorsque la taille maximale de la fenetre tend vers l'infini. 

5 De plus / la presente invention permet d' analyser les fluc- 
tuations instantanees et aleatoires du debit, ce qui peut 
etre utile pour estimer la quality de service offerte ci une 
connexion* Elle s'adapte aussi bien pour des simulations 
efficaces de la dynamique d'une session TCP operant sous 
10 controle, de bout-en-bout, sur un reseau de grande dimension o 

Ci-apres, il est donne une representation max-plus generale 
du module de base qu' utilise la presente invention* Ensuite, 
les services deterministes avec une Evolution deterministe 
puis markovienne de la taille de la fenetre sont decrits . Le 
Demandeur a montre que pour ces modules deterministes, le 
debit ne depend que du temps aller-retour "round trip time" 
ou RTT et de la bande passante (ou "bandwidth"). Des varian- 
tes et extensions du module deterministe ont ete aussi mises 
en oeuvre, notamment s 

- le cas des pertes aleatoires en plus des pertes de conges 
tion; 

- un suivi detailie de la taille du buffer et des pertes de 
congestion; 

- le cas de services aleatoires; 

- le cas en presence de time-outs* 

Les extensions du modele deterministe conduisent toutes a des 
30 formules analytiques ou ci de nouveaux principes de simula- 
tions bases sur le calcul du produit d'un grand nombre de 
matrices. En particulier, le Demandeur a montre que le cout 
de la simulation par cette approche de la transmission de n 
paquets sur K routeurs est avantageusement en 2n(KW max ) 2 , oil 
35 w max est la faille maximum de la fenetre. 

L'"algdbre" max-plus scalaire est un semi-anneau sur la 
droite reelle ou 1'on remplace 1' addition par max (notee ©) 
et la multiplication par plus (notee ®). La loi ® est 



15 



20 



25 
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distributive par rapport a la loi ©, ce qui permet d'etendre 
les concepts usuels de 1' algebre lineaire & ce cadre, et en 
particulier la theorie des matrices. Ce semi-anneau est note 
(R maxf © f 0) , ou est la droite reelle completee avec moins 

5 l'infini, qui est 1' element neutre de ®. Par la suite, 
1' ensemble des matrices carrees de dimension d dans cette 
algebre est note (1^^,©,®), otl les deux operations © et ®, 
lorsqu'elles sont appliquees a des matrices, sont liees par 
la relation El de 1' annexe I. 

10 

Pour plus de details sur cette algebre, qui est egalement 
utilisee pour des garanties QoS dans les reseaux, on pourra 
se referer a [2] ou [6] dans les references donnees & 
1' annexe II. 

15 

Dans un reseau PAPS, de type "premier arrive, premier servi", 
avec K files en tandem, le n-ieme client arrivant & la 
station i regoit un service a ± (n) . Dans le contexte de TCP, 
ce reseau modelise une seule source envoyant des paquets vers 

20 une seule destination, & tr avers un chemin compost de K 
routeurs. La variable a ± (n) est le retard aieatoire cause par 
le trafic transversal (les autres utilisateurs ) au niveau du 
routeur i sur le n-ieme paquet. Ce retard n'inclut pas les 
delais d'attente, mais seulement le ralentissement de la 

25 vitesse du serveur du h la presence du trafic transverse. Le 
delai de propagation du n-ieme paquet entre les routeurs i h 
j est note ci-apres d i,j<n). 

Le debit du flux d' entree est controie par une fenetre 
30 dynamique W(n), dont la taille est egale au nombre total de 
paquets envoyes par la source a un instant donne et n'ayant 
pas atteint la destination (ou plus precisement les paquets 
qui n'ont pas encore ete "acquittes" ) . 

35 La taille de la fenetre a une evolution generale definie par 
la recurrence E2 donnee en annexe. 
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Dans les Equations E2, ACK(n) est le signal de contr61e de 
flux/congestion donnant 1 • information sur l'etat du systeme 
au temps n et ou f est une fonction precisee ci-apres. 

Par exemple, ACK(n)=l si aucune congestion (ou aucune perte 
de paquet) n'est observee par le n-ieme paquet, sinon 
ACK(n)=0. Par la suite et dans certains cas, f est supposee 
dependre aussi de Ws, qui est, pour les modeles TCP, le seuil 
qui separe la phase de croissance exponentielle (dite "slow- 
start phase") de la phase de croissance lineaire (dite 
"congestion-avoidance phase"). 



La relation de recurrence (donnee en annexe par E3 ) definit 
une taille de reference de la fenetre. La taille effective 
15 est alors d«§finie comme la partie entiere de la taille de 
reference, selon 1' Equation E3 de 1' annexe. 

La taille maximale de la fenetre est supposee finie et notee 
w* (voir 1' equation E4 de 1' annexe). 



De plus, devolution de la taille de la fenetre peut etre 
decomposee en deux phases qui dependent de ACK(n) de la facon 
suivante: une phase de croissance definie par E51 et une 
phase de decroissance definie par E52. 

Dans le modele d^crit ici, ACK(n) (valant 0 ou 1) est le 
signal d ' acquittement du n-ieme paquet qui detecte un etat de 
congestion ou une perte de paquet. Les exemples usuels de la 
politique ideale de 1' evolution de la fenetre sont donn#s par 
les equations E61 a E67 de 1' annexe, dans lesquelles le reel 
a est compris entre 0 et 1, strictement. Les relations E61 et 
E62 sont donnees respect ivement pour des phases exponentielle 
et lineaire. 



35 Dans les exemples suivants, la valeur de a est fixee a 1/2. 

La file a 1' entree est consideree ci-apres comme saturee 
(meme si le cas non-sature peut etre integre facilement, 
comme on le verra. plus loin). Le reseau se comporte alors 
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coirime un reseau ferme et son debit donne le taux maximum 
auquel la source peut envoyer des paquets tout, en gardant un 
buffer d' entree stable. On pourra se ref6rer a 1' article [5] 
donne en annexe des references pour plus de details sur ce 
cas . 

On se refere a la figure 1 montrant K files en tandem, avec 
controle de flux. x ± (n) est la date a laquelle le n-ieme 
paquet arrivant au routeur i commence son service sur ce 
routeur . 

y ± (n) est la date a laquelle le n-ieme client quitte le 
routeur i. Par les Formules E7 de 1' annexe, on en deduit le 
vecteur Z(n), dit "vecteur des dateurs" par la suite. 

Les variables Mi, i appartenant a {1, w*}, sont des 

matrices donnees de (F^4^). e est la matrice de (I^ax) dont 
tous les elements sont egaux h moins ljinfini. Ci-dessous, on 
note (M 1 |M 2 | . . . 1^*) la matrice de (I^L'^ > d^finie par blocs 
de taille k x k. Le Demandeur a montre que tous les blocs 
sont egaux & la matrice € de (I^*) ' sauf pour la premiere 
ligne de blocs qui est egale & M X ,M 2 , . . . ,M W * . Le Demandeur a 
montr6 aussi que si initialement le systeme est vide, 
1' evolution du vecteur des dateurs Z(n) est donnee par la 
recurrence max -plus lineaire E8. 

Dans la formule E8, D est la matrice de dimension Kw* dont 
tous les elements sont egaux h moins l'infini sauf ceux de la 
forme D K+i ± (avec i appartenant a {1 , . . . , k( w*-l ) } ) qui sont 
tous egaux h 0. 

On peut definir l'exposant de Lyapunov de la station k par la 
limite de 1' expression y k (n)/n quand n tend vers plus l'infi- 
ni. A cause de la monotonie, il est clair que cette limite 
est independante de k. Cette propriete est plus generalement 
vraie sous l'hypothese d' irreductibilite definie dans la 
reference [8]. 
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Au niveau de la representation decrite ici, le controle 
d'erreur n'est pas pris en compte et aucune difference n'est 
marquee entre les paquets originaux et les paquets retrans- 
mis. En particulier, il n'est pas encore etabli de distinc- 
5 tion entre le "send rate", throughput" ou "goodput" au sens 
de la reference [13]. 

L' equation E8 est, dans cette realisation, la base du schema 
de simulation algebrique auquel il etait fait allusion ci- 
10 avant. Comme les matrices A^^n) sont de dimension Kw*, et 
comme seuls les produits matrice-vecteur sont necessaires, on 
peut simuler la transmission contrdlee de n paquets a travers 
le reseau en 2n(Kw*) 2 operations sur un seul processeur. 

15 La figure 2A represente un exemple d' evolution pas-a-pas 
explicite des dateurs et de la taille de la fenetre (K=5 avec 
io lf ....o 5 ) = (1, 1, 2, 1, 1) et w° = 4). revolution de la 
taille de la fenetre est deterministe (TCP Tahoe sans la 
phase de "slow start", soit (w 1# w 2 , w 3 , ) = (1, 2, 2, 3, 

20 3, 3, 4, 1, 2, 2, ...)). 

Par exemple, le paquet # 5 correspond a la taille de la 
fenetre 3,'ce qui est note 5(3) sur la figure 2A. Ainsi, 
juste apres la transmission du paquet # 5, trois paquets dans 
25 le reseau n'ont pas encore ete "acquittes". Comme le montre 
la figure 2A, les retards imposes aux paquets (voir par 
exemple le paquet # 6) dans les routeurs internes peuvent 
§tre assez compliques. 



30 



Pour la simplicity de la presentation, on considere d'abord 
le cas oii tous les delais de propagation d if j(n) valent 0, et 
ensuite une adaptation pour couvrir le cas de delais de 
propagation non nuls . 



35 



Dans une premiere approximation, les services o ± (n) sont tous 
deterministes et constants (a i ) . La congestion est normalement 
detectee par la politique de la formule E9. 
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Cette detection peut §tre interpretee comme donnant le taux 
de service 1/a* du goulot d' etranglement du reseau et S, un 
temps d' aller-retour RTT. On detecte done une congestion 
quand le taux d' Emission moyen (w n /RTT) atteint le taux du 
5 goulot d' etranglement . Pour le modele a fenetre statique, w* 
est la taille de la fenetre optimale donnee par le produit : 
bande passante x delai (reference [9])- 

Le Demandeur a montre que, dans ces conditions, w* est donne 
10 par la formule E10 et que w n devient periodique avec une 
periode T qui peut etre decomposee en t x +t 2 + . • • +t w * , ou ti est 
le':- nombre d' occurrences de w n = i pendant une periode. 

En considerant d'abord le modele TCP Tahoe sans la phase 
15 exponentielle, sous la politique E9 , il peut etre deduit les 
formules Ell de 1' annexe* La matrice H (formule E12) est 
irreductible et sa valeur propre est egale h nS+a* . Le 
Demandeur a aussi montre que si tous les delais sont deter- 
ministes, cette valeur propre est donnee par la formule E13. 

20 

Par consequent, le debit depend des delais (a 1# ...a k ) seule 
ment & travers S et o* et ce debit est donne par la formule 
El 4. Ainsi, le debit asymptotique devient l/(2a*) quand w* 
tend vers plus 1'infini. 

25 

Sur la figure 2 (ou l'on suppose que w* = 6 ) Involution pas- 
a-pas de la date d' entree y 0 (n) et de la date de depart y k (n) 
(du paquet # n) sont representees. Les proprietes algebriques 
ci-avant s ' interpretent comme suit. 

30 

Avant la detection de congestion, les paquet s envoyes se 
comportent comme s'il n'y avait aucune interaction entre eux, 
sauf pour le couple de paquet s envoyes en meme temps lorsque 
la taille de la fenetre croit d'une unite. Pour ces derniers 
35 paquets, le second paquet quitte tou jours la station K avec 
un d£lai de a* par rapport au premier. On peut done lire 
directement la valeur propre 5S+6a* sur Involution graphique. 



SDOCID: <WO 0165772A1 1 > 



WO 01/65772 PCT/FR01/00579 

17 

Le d6bit de saturation pour le TCP Tahoe avec la phase 
exponentielle est donne par la formule E15 et le debit 
asymptotique (w* infini) est donne par E16. 

5 La figure 3 represente 1' allure du rapport n/y 4 (n) (traits 
pleins) et de Wn (traits point illes) dans les conditions 
suivantes : TCP Tahoe sans la phase exponentielle avec quatre 
files en tandem (0^=3, 2; o 2 =4,61; o 3 =2,7; o 4 =4,61 et w*=4 ) . Le 
d^bit de saturation est egal a 0.140084 en utilisant la 
10 formule E14. 

Une evolution periodique deterministe de TCP Reno a ete 
consideree dans la reference [11] pour obtenir une valeur 
heuristique du d6bit. La representation max-plus ci-dessus 
15 conduit a une nouvelle Formule qui raffine celle de [11]. 

Si tous les d^lais sont deterministes , le d6bit de saturation 
ne depend que de S et de o* et il est donn6 par les formules 
E17 de 1' annexe. 



20 



25 



En se r6ferant a la figure 4, les debits asymptotiques 
(cas 1, cas 2, cas 3) s'obtiennent de maniere assez intuitive 
a partir de 1 ' approximation fluide de 1' evolution de la 
taille de la fenetre. 



On appelle d 0 =l/a* le d^bit correspondant a w n = w (cas 0) ; 
lorsque w n augmente lin^airement a partir de 1, le volume du 
debit, qui est proportionnel a l'integrale de W(t) sur une 
periode, est bien 1/2. d Q (cas 1); lorsque w n augmente 
30 lin^airement a partir de w*/2, le volume du d6bit diminue d'un 
facteur de 3/4 (cas 2). 

La formule bien connue du d6bit d'une connexion TCP en 
fonction du taux de perte P per te et du tem P s aller-retour RTT 
35 est de la forme donnee par la reference [11] et retranscrite 
par % 

dSbit = C C /(RTT. (P p erte) 1/2 ) 

ou c c est une constante r6elle. 
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Pour le modele dEterministe ci-avant , avec RTT^S, les valeurs 
de p perte selon les differents cas (CI, C2, C3) s'expr invent 
selon les formules E18. 

Par consequent, pour les grandes valeurs de w* (ou les petites 
valeurs de P perte )' la fondle asymptotique coincide bien avec 
celle de la reference [ll], 

Le Demandeur a verifie que tous les rEsultats ci-dessus sont 
vrais pour des dElais de propagation d if j constants a condi- 
tion de remplacer la valeur de S par celle de la formule El 9 
de 1 ' annexe . 

Pour les modeles deterministes, le d<§bit obtenu peut 6tre 

15 compare a celui simule par un simulateur (par exeraple de type 
NS) en choisissant une taille de paquet arbitraire et en 
prenant une vitesse de service du routeur i correspondant & 
a ± . Les ecarts entre les debits obtenus par la simulation NS 
et par les formules ci-avant ne peuvent provenir ici que de 

20 differences sur les mEcanismes de detection de la perte par 
congestion, ou du fait que la partie entiere de W(n) dans f 
(comme dans la reference [9]) est seule considErEe ici, 
tandis que NS utilise W(n). En effet, pour tous les modeles 
deterministes avec la meme Evolution periodique de W(n) , les 

25 Evolutions sont exactement les memes. Par exemple, un 
protocole TCP simulE sur NS avec une source de type ftp 
donne, pour K=10, une taille d'un paquet egale & 1250 (40 
pour les accuses de reception), et une taille des buffers 
Egale k 2. Tous les delais d ± ^ sont egaux ^ 0.1ms sauf d K0 

30 qui est Egal a 1ms. La vitesse de service vaut 
(10,5,4,2,5,4, 5,5,4,5,5)Mb/s pour les liens 0-1, . . ,9-10, 10-0. 
A t=100s, le simulateur NS donne 152.27 paquets/s. Pour cet 
exemple, S=25,5 ms et a*=5 ms; done, par la formule E15 on 
aurait 134 paquets/s. Cependant, w* est en fait Egal & 7 dans 

35 la simulation NS au lieu de 6 dans le modele dEcrit ici. Un 
RTT est nEcessaire pour dEtecter une temporisation dite 
"triple-acks" sur ce module NS. En Evaluant le dEbit avec la 
formule E15 et w*=7, on a 152.55 packets/s. 
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Enfin, pour tous les modeles deterministes , la suite {w n } est 
dyterministe, et on peut done calculer Involution du dybit 
d' Emission de maniere exacte & partir des produits de 
matrices max-plus correspondants • 

5 

On d6crit ci-apres un modele de Markov avec services d6ter- 
rainisteso 

Les services sont consid6r6s tous comme dyterministes • 
10 Maintenant, W(n) est suppos^e yvoluer ind^pendamment des 
autres yiyments du r^seau par la matrice de transition 
markovienne P donn6e par les formules suivantes s 

- pour tout n entier, 

15 

o si w ( n ) <w max , w(n+l) = 

w(n) avec probability p Q 
h(w(n)) avec probability p+ 
g(w(n)) avec probability p. 

20 

. si W ( n ) =W max , W(n+l)=g(w(n) ) avec probability 1, 

oil h(W(n) )>W(n) , g(W(n))<W(n) et oil W max est arbitraire, 

25 Si tous les services sont dyterministes , pour toute yvolution 
markovienne de W(n) avec la structure ci-dessus, le dybit ne 
dypend des delais qu'& travers S et a*. De plus, l'exposant 
de Lyapunov y ( 1 ' inverse du dybit) est de la forme donnye par 
la formule E20 dans laquelle n est la probability station- 

30 naire de 1' ensemble A. 

Dans le cas ou les deux transitions sont possibles, w* peut 
etre dyfini par la politique E9 et W max =w*, avec la politique 
E21 de 1' annexe. Ce module peut etre compary h celui de la 
35 ryfyrence [12], ou une probability de perte globale est 
utilisye pour capturer en meme temps les time-outs (TO) dus 
aux pertes des paquets et les triple-duplications des ACKs 
(TD) dus h la congestion. Dans le module dycrit ici, ces deux 
mycanismes sont au contraire dycrits s6par6ment; les pertes 



INSOOOO: <WO 016S772A1 I > 



WO 01/65772 PCT/FR01/00579 

20 

de paquets qui generent les TO constituent une suite i.i.d. 
(independante et identiquement distribute ), independante de 
tous les autres Elements du reseau et sont capturees par le 
parametre p_; les pertes dues a la congestion sont capturees 
5 par le paramfetre w*. 

On montre que pour TCP Tahoe avec une evolution markovienne 
comme ci-dessus (c'est-a-dire avec E21 et g(W(n))=l), le 
debit est donne par les formules E22. 

Le Demandeur a verifit que p. = 0 correspond au cas de TCP 
Tahoe deterministe (sans la phase exponentielle) • Par 
ailleurs, si les pertes de paquets sont liees a la congestion 
d'une fagon qui est en moyenne similaire a celles du cas du 

15 TCP Tahoe deterministe , le parametre C Q vaut 60%. Done, la 
degradation des performances peut etre significative (ap- 
proximativement 15%) lorsque qu'on passe du modele determi- 
niste (ou C o =0,71) au modele markovien. L ' impact des TO sur 
les performances a deja ete remarque par exemple dans les 

20 references [13,10]. L ' influence preponderante de p- sur la 
probability de perte globale peut 6tre quantif iee h partir du 
modele analytique ci-avant comme suit. 

La probability de perte globale de ce modele est donnee par : 

25 

/J<w*)+p_(l-ji(w*))=ji(l) 
ou jj(w*) est la perte due a la congestion et p_(l-jj(w*) est la 
perte due aux TO. 

30 Pour fixe, la Figure 5 montre la decroissance du debit 

en p_ (a*=l et S=w*-1 ) . 

On voit Ik dans quelle mesure 1' impact des pertes par TO est 
preponderant par rapport aux pertes dues & la congestion. 

35 

En appliquant les modeles ci-avant aux memoires tampons (ou 
"buffers") d'un reseau, des resultats similaires sont obtenus 
pour des modeles du type TCP Reno et aussi pour des modeles 
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bas^s sur une evolution markovienne de la taille de la 
fenetre avec trois transitions (ou plus). 

Par exemple, soit K=4 avec 0^3,2} a 2 =4,61; o 3 =2,7? a 4 =4,61, 
5 et w *=4. La Figure 6 montre Involution de n/y 4 (n) et w n pour 
TCP Reno markovien avec (p+/Po'P~) = (0.8,0.1,0.1). 

La Figure 7 montre les debits de TCP Reno et TCP Tahoe dans 
les cas d^terministe et markovien ( (p + ,p 0 1 P- ) = ( 0 • 8 * 0 o 1 , 0 . 1 ) ) . 
10 Dans le cas d^terministe, les debits sont 6gaux & 0-140084 
(Tahoe) et 0.172166 (Reno). Les debits du module d^terministe 
semblent sup6rieurs h ceux des modeles markoviens. 

Pour des buffers finis, des files en tandem avec des services 
15 d^terministes et avec une taille du buffer fixe b pour toutes 
les stations sont d'abord consid6r6es. Dans ce cas, d'aprfes 
la formule de Little, il est naturel de definir la detection 
de la congestion par la formule E23o 

20 Ce modele conduit encore & un comportement d^terministe et 
periodique de W n et & une analyse similaire a celle d^crite 
ci-avant. On remargue que le modele de la politique E9 est un 
cas particulier de E23 avec la taille du buffer b=lo 

25 Pour des services al£atoires, des files en tandem avec des 
services al^atoires et avec une taille du buffer fixe b pour 
toutes les stations sont maintenant consid6r€es. Les supports 
des temps de service sont supposes born6s , 

30 On considfere alors les relations E24 exprimees en annexe, 

W(n) a une Evolution d6terministe lorsque W(n) appartient a 
{l,...,b} et une Evolution al^atoire lorsqu'il appartient & 

{b+1, ,b+k}. L' interpretation de ce module est la suivantes 

35 pour une Evolution de la taille de la fenetre du type TCP, la 
date de d6but de la congestion est le moment ou 

w(n) = S(n)/a*(n) 
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Ensuite, une perte de paquet a lieu lorsque le nombre de 
paquets dans le systeme d^passe la taille du buffer , ce qui 
correspond sensiblement a l'6v£nement W(n) = b + r(n). 

Une autre interpretation est la suivante: b est une taille 
arbitraire (par exemple, la moitie de la taille r^elle du 
buffer) telle que si le nombre des paquets en attente d6passe 
ce seuil, la congestion est d6tect6e. 

Sous la politique E24, si les vecteurs de temps de service 
o if ...a k sont i.i.d. en n, le debit de TCP Tahoe sans phase 
exponentielle est donne par 1' expression E25 ou 0 est le 
vecteur (0, ,0)* de R k ' 1 - 

Les mecanismes de time-out peuvent etre pris en compte, par 
exemple, par la condition S(n) > TO ou par la condition: 

Yk< n ) ~ YK< n -n-i> > TO - 

Par exemple, des que Z(n) est markovien et que le support D 

des valeurs absolues {y K ( n )-YK( n - w n-i ) > est fini, on retrouve 
la meme structure regenerative comme ci-dessus. La loi d'un 
cycle de regeneration Tl peut etre explicitement calcul^e par 
la recurrence E2 6 de 1' annexe. 

Cette formule de recurrence est valide pour n<W max . Ainsi, on 
peut obtenir une formule pour le debit en utilisant le 
th^oreme ergodique pour les processus reg^neratifs comme ci- 
dessus. 

L'homme du metier peut obtenir une formule similaire pour TCP 
Tahoe avec phase exponentielle ou TCP Reno, et aussi pour 
diverses extensions du modele ci-dessus en incluant des 
pertes de paquets comme suivant la formule E21 ci-avant. 

Dans un autre exemple, ou 1'on considere guatre files en 
tandem avec une taille du buffer b=50 et avec une distribu- 
tion multinomiale independante pour les services, la figure 
8 montre 1' evolution de n/y 4 (n) et W(n)pour le modele TCP 
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Tahoe avec la phase exponent ie lie et avec b = 50. Les o ± (n) 
sont i.iodo et mutuel lement ind<§pendantes , ^ valeurs dans 
{1,5,10} avec des probability respectives de 0.3, 0o4, 0.4o 

La figure 9 montre une comparaison des debits de deux TCP 
Tahoe avec la phase exponentielles b = 10 et K=10. Dans le 
premier cas, a 3 (n) = l et {o^n) (avec j*3)> est une suite 
i.i.d., h valeurs dans {1, 10, 20} avec les probability 
respectives 0*3, 0.2, 0.5. Dans le second cas o 8 (n), est 
constant <§gal h 1 et les autres services sont comme ci- 
dessus. 



Le calcul des debits au moyen de la formule E25 montre que la 
permutation des caract6ristiques de deux routeurs peut 
15 changer la valeur du d6bit- Done, d<§ja dans ce cas, le module 
ne peut 6tre r€duit h un module avec un seul serveur, celui 
du goulot d'6tranglement . ' II convient de noter aussi que la 
connaissance des valeurs moyennes est insuffisante pour 
pr^dire la valeur du d6bit moyen. 



Le module d6crit ci-avant trouve de nombreuses applications , 
notamment s 



- relativement a la loi des debits instantan^s qui, dans le 
25 cadre de la formule E20, est 6gale h 1' expression E27 de 
1' annexe ; e'est une valeur importante qui peut d^finir un 
indicateur naturel de QoS en complement de la valeur moyenne; 



30 - la loi des d^lais bout-en-bout (formule E28) ; 

- la loi du temps T n6cessaire a la transmission d'un fichier 
de taille F; en premiere approximation, cette loi est donn6e 
par la relations P(T>t) = P(y K (F)>t) ,mais une formule plus 
35 precise peut etre donn^e en prenant en compte des retrans- 
missions de paquets perdus . 

Dans ce cadre, les problfemes suivants peuvent Stre trait6s au 
moins dans le cadre de la simulation alg^briques 
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- les modeles ouverts (ou la source n'est pas saturee: HTTP 
au lieu de FTP) ou le processus des arrivees est d£crit par 
ses caracteristiques statistiques; 

5 - les connexions multiples et les interactions entre plu- 
sieurs utilisateurs ; 

- les connexions multipoints a travers un rgseau avec une 
structure arborescente au lieu d ' une structure lin^aire de 

10 routeurs en serie (voir la reference [7] pour le cas a 
fenetre constante). 

Quelques precisions sont donnees ci-apres sur la structure 
d'un simulateur algebrique qui permettrait de pr§voir les 
15 performances de TCP ou d'autres mecanismes de controle de 
flux issus de TCP dans des reseaux IP existant d£ja ou en 
cour s de conception • 

Le principe de base de ce simulateur est 1' analyse d'une 
20 connexion contr616e. La simulation de cette connexion 
commence par 1 ' acquisition des donnees concernant le r^seau: 

- nombre des routeurs traverses par la connexion/ carac- 
teristiques de chacun des routeurs: capacity en Mb/s f 

25 taille des buffers bi, previsions, etc; 

delais de propagation sur chacun des liens (d ± i+1 ); 

caract£ristiques statistiques des flux transverses F ± 
30 lorsque ces dernieres sont connues (comme par exemple 

dans le cas des routeurs du backbone Internet); 

caracteristiques statistiques du traf ic of f ert par la 
connexion controlee; 



35 



caracteristiques statistiques des traf ics of ferts par 
les autres flux en competition avec la connexion contrd- 
lee sur le routeur d'acces. 
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Ces caract^ristiques apparaissent sur la figure 10. 

On construit alors en fonction de ces donnees et des carac- 
t^ristiques de la version de TCP choisie (ou de 1 ' algorithme 
5 de controle de flux retenu) un simulateur alg^brique qui 
consiste en la construction pas a pas des matrices ci- 
avant et le calcul par recurrence des dateurs du systfeme 
comme indiqu6 • 

10 On d^duit de 1' analyse math^matique qui a 6t6 pr€sent6e ci- 
avant que l'on peut calculer a partir de la simulation les 
principales caract<§ristiques de l'effet de TCP sur la 
connexion contr616e (d£bit moyen, d6bit moyenn6 sur une 
p^riode de temps donnee, fluctuations du d6bit instantan^ 

15 etc). 

Pour completer cette description du simulateur, quelques 
indications sur la maniere de prendre en compte certains 
m^canismes fins sont donnees ci-aprfes. 

20 

Pour un routeur i, soit bi la taille de son buffer par flux; 
la detection du d^passement de cette capacity se fait par la 
condition suivante sur les dateurs; 

x^n+b^syifn) 

25 

Comme deja indiqu6, les m6canismes precis des time-outs 
peuvent etre pris en compte par la condition ; 

yK< n > - YK( n - w n-l) > T0 

30 Cette formulation vaut aussi lorsqu'il y a Evolution de TO 
lui-meme (dans certaines versions, la variable TO est mise h 
jour h chaque timeout: TO est double en cas de timeout, 
jusqu'& ce qu'on atteigne 64 fois la valeur d' initialisa- 
tion). La pr^sente invention peut s'appliquer & 1' estimation 

35 des variations et aux adaptations de la fenetre dans le cadre 
max-plus lin^aire. Les variations et adaptations de la 
variable TO sont de nature semblable et sont done aussi 
repr^sentables dans ce cadre. 
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Dans certaines versions de TCP, il est recommande de proceder 
au regroupement des accuses de reception. On attend par 
exemple que trois paquets soient arrives a la destination 
pour renvoyer un seul accuse pour les trois paquets. Ceci 
5 peut aussi etre represents au moyen d'un chalnage de deux 
mecanismes de fenetre variable, et rentre assez naturellement 
dans le cadre d'une description max-plus lineaire. 

Un mecanisme de scheduling dans un routeur (par exemple 
10 premier arrive premier sorti "FIFO" (ou Weighed Fair 
Queueing) avec eventuellement des dif f erenciations de 
services) peut etre represents dans cette simulation par la 
prise en compte de 1' influence de ce mecanisme sur les durSes 
de service a ± (n) ; une fois calculees en fonction du trafic 
15 transverse, du trafic de la connexion controlee et du 
mecanisme de scheduling, ces durees aleatoires peuvent etre 
injectees dans le simulateur. 

Pour des infrastructures de reseaux d'acces en cours de 
2 0 conception, on ne peut evidemment pas mesurer les caracte- 
ristiques statistiques des trafics transverses sur les 
routeurs d'acces. En supposant qu'on dispose d'une descrip- 
tion statistique du trafic offert prevu sur les routeurs 
d'acces de ce reseau (par exemple N trafics HTTP), on peut 
25 alors utiliser le simulateur pour calculer le trafic trans- 
verse inconnu coirane un point fixe du systeme: soit t la loi 
du trafic de la connexion controlee observe au niveau du 
routeur d'acces; t doit etre tel que si on prend coitime trafic 
offert sur la connexion controlee celui d'une session HTTP, 
30 - le trafic de la connexion controlee observe sur le 
routeur d'acces a pour loi t; et 

- le trafic transverse observe en ce point est la somme de 
N trafics de loi t. 

35 Ainsi, dans les reseaux satures ou non, le mecanisme de 
controle en feedback adaptatif de TCP est un feedback 
lineaire dans l'algebre max-plus. Ceci conduit a une repre- 
sentation simple de l'effet de ce protocole sur un reseau 
quelconque qui admet lui-meme une representation max-plus en 
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1' absence du controle, comme c'est le cas par exemple pour 
les files en tandem ou les r<§seaux fork- join que l'on trouve 
dans les arbres multipoints. 

5 Les formules ci-avant confirment que dans ce cas, le d6bit ne 
depend que du RTT et du taux du goulot d' Strang lement- De 
nouvelles formules ont aussi 6t6 obtenues lorsque les 
services sont al<§atoireSo Les aleas repr^sentent ici l'effet 
du reste du trafic sur la connexion contr616e- Dans ce cas, 
10 le dSbit ne peut pas etre obtenu seulement h partir de 
considerations en moyenne, et que l'ordre et le comportement 
statistique fin des routeurs ne peuvent etre ignores. 

L' ensemble des modeles entrant dans ce cadre est tres riches 
15 - on peut en effet choisir des services deterministes ou 
aieatoires, un controle de flux base sur la perte ou conges 
tion; 

- les pertes peuvent provenir de la congestion ou des time- 
outs, ou etre aieatoire, ou encore etre une combinaison de 

20 ces trois possibility y 

- une version de Reno ou Tahoe peut etre choisie, avec ou 
sans la phase exponentielle etc. 

La presente invention trouve des applications interessantes 
25 dans 1' analyse de ces quelques combinaisons , puisque toutes 
les combinaisons peuvent etre en principe analysees dans ce 
cadre o 

Plus generalement , 1' invention definit un cadre g6n6rique 
30 pour la simulation des protocoles du type TCP sur des rSseaux 
qui peuvent etre grands . La simulation est bas€e sur un 
traitement simple qui exploit e ia linear ite et qui a une 
complexity maltrisee.. 

35 Bien entendu, 1' invention ne se limite pas a la forme de 
realisation decrite ci-avant & titre d' exemple; elle s'6tend 
h d'autres variantes qui restent neanmoins definies dans le 
cadre des revendications ci-aprfes. 
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E4 : {AGB^^AikGBkj^maxiAik.Bv), (A ® B)^ — 0 A« 9 B ki = max A* + B w . 

i<*<* 

£2 : ^(0) = l, W(n + l) = /(W(n)MCJf(n)), 



E 54 : {n:^CJT(n) = l>, 0</(W(n).D-W(»)<l. 

£ ^2 : {n : ACK (n) = 0} ' 1 ~ '^WOT* w(n) " 

TCF^Tbhoe: 

. /(W(n),l,W.(n)) = W(n) + 1; W(n)<W.(n), 

„ ^ = W[rx) + — ; W(»)>W.(n), 

, W,(n + 1) : = W.(n); 

; /(W(n),0,W.(n)) = i ; 

£6 5 ; ^-(n + 1) : = laW(n)j . 



Td* Reno 

BC6 : Wn),!) = wr(n) + -i-, 

IDfi 

Efi> : /(W(n),0) = l<*^(n)J. 

£ ^ : V « > o, ff °< n ) = 0 ; = + »«(») 

yo(n) = yjr(n — w»-t) ® d Kfi {n - iu„-i), 

W(n) = [itt-i (") ® *-i © V« (r» - 1)] ® » = 1 » • • • • K • 

= (yi(n).w(»»).-»*( n )) € 

S(n) = (K(n),r(n- l),,.,y(n - u>' + 1))' € 



i/VO 0165772A1 I > 



WO 01/65772 



PCT/FR01/00579 



29 



£ Q ; Z{0) = (0, ... 0)«, Z(n) = (n) 8 Z(n - 1), Vn > 1, 

= (M(n) ® AfJ(n)|£|.-|£) ©D, ^ a (n) = (M(«)|Mi(n)|C|..|f ) © D. 
(M(n)|f|..|f|Afi.(n))©D. 

Ki/W : («/(»))« = E < d »-» t(n) + "» (n)) + dK ° (n " 0 [ i = J 3 J "°° P < ^ 
£^0." ti;* =min{n:a'n>S}= [^-J+l<K + l. 



E 4d: 



{wi,u>a,...} = {1,2,2,3,3,3,... - - 1, it;', 1,2, 2,-). 



Vi e {l, ... »• -i}.** = <,-.*-• -l; r = + 1 

Vi > 1, 

/AfffiAf' ATM* M*~ 1 M' £ --- \ 



A/'" 1 Af' 
: f 
Af» 
AT 



* • * 

e at 

£ 



V»>1 , Vn 6 {1, 1). 



/ M n £ ••- £ Af' 
: : £ 

M 

la £ £ 
£ l 4 £-£ 



M n ~ l M' £-\ 

Af' 
£ 
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H « M ® (M n ffl Af') ® • • • ® (hi ffi M*) 

E >13 A = (w*-l)S + a\ 

IT i; , - l) + 2 

fc. A 1 * '• 7 ~ 2(W-l)S + <r»" 

P ,C. T -,_l ^C-l)-|f|(lV|-3) 



E M£ or 



• 1 _i ^'-i)-L^jg^j-i) ^2 

2 (»--[Hlj)S + < r. • 



11T* - . _^ _i 3 1 



C2 : i)-iV 



-i)-m(L^J-i)+2- ' 



t . A = {S — (* — l)a*fk = 1, .., i»* — 1} u {a*}, 

ft - {»Jr(" + 1) - yjf (n)}»ew 

■ MWW) = ^{n) + -i . ^ = 0. 
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£2.3: ACK(n) = l, if »„<6+ JL, 0, ifti»n>6+— . 



m 1 = " ax 'i^' i|n) 



S = sup{5(n)} ; a' = sup{a»(„)} ; r(n) = ; «" = 6 + nup{r(n)} 

ACK{n) = 

It vh, < b + r(n), 
°» ««»>» + r( n ). 



P(T, > n, 2„ = x) = P(n: =1 {yx(*) - yjr (* - ) < TO, Z n = *) 

5T ^("2=1 {l//f (*) - yjc(* - it;,.!) < TO, = r', = z) 



CCA 

= 23 P(Z„ = z, yjf(n) - y/f («/„_,) < TO | Z„_, = z')P(r, > n - 1, = 
*. Bmn— i» (&,*■(„ + i) _ 9te(n)) -i < « _ , 



; lim ir" 



Km — "» & J» (yir(n) , y K { n - «»„_,) < x) 
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Re vend ic at ions 

1. Dispositif d'aide a la surveillance et/ou de simulation 
d'un systeme complexe, notamment d'un reseau de communica- 

5 tion, caracterise en ce qu'il comprend: 

- une memoire ( ROM , RAM) pour stocker de premieres donnees 
representatives de parametres du reseau, ainsi que pour 
recevoir au moins de secondes donnees representatives 
d'evenements dans le reseau, une portion de ladite memoire 

10 etant reserv£e au stockage de donnees sous forme matricielle, 

- un module de calcul (CAL), apte a effectuer sur au moins 
deux matrices de structure dynamiquement variable, une 
operation formant produit selon l'algebre dite MAX-PLUS, 

- un module de modelisation (MOD) pour construire au moins 
15 une premiere matrice (2(n)) et une matrice courante (A(n)) 

respectivement en fonction des premieres donnees et des 
secondes donnees , selon un modele choisi, et 

- un module pilote (PIL) pour appliquer repetitivement la 
premiere matrice et la matrice courante au module de calcul, 

20 la matrice produit obtenue <Z(n+l)) devenant une nouvelle 
premiere matrice. 

2. Dispositif selon la revendication 1, caracterise en ce.que 
les premieres donnees comprennent des donnees relatives h la 

25 topologie du reseau. 

3. Dispositif selon la revendication 2, caracterise en ce que 
les premieres donnees comprennent le nombre (K) de routeurs 
du reseau traverses par une connexion a surveiller ou a 

30 simuler et des proprietes desdits routeurs. 

4. Dispositif selon la revendication 3, caracterise en ce que 
les premieres donnees comprennent des tailles memoire des 
routeurs (b K ) . 

35 

5. Dispositif selon l'une des revendications 2^4, caracte- 
rise en ce que les premieres donnees comprennent des proprie- 
tes statistiques des trafics offerts dans le reseau. 
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6. Dispositif selon l'une des revendications precedentes, 
caracterise en ce que le module de modelisation comprend : 

- un sous-module de modelisation statique (ST), pour cons- 
truire la premiere matrice en fonction des premieres donnees, 

5 et 

- un sous-module de modelisation dynamique (DYN), pour 
construire au moins une matrice courante en fonction des 
secondes donnees. 

10 7. Dispositif selon l'une des revendications precedentes, 
caracterise en ce que les secondes donnees comprennent des 
informations relatives a des pertes (p__,p+) dans le reseau. 

8. Dispositif selon l'une des revendications precedentes, 
15 caracterise en ce que les secondes donnees comprennent des 
informations relatives a des flux transverses dans le reseau, 
par rapport a une connexion contrdlee a surveiller ou h 
simuler. 

20 9. Dispositif selon l'une des revendications 7 et 8, carac- 
terise en ce que les secondes donnees comprennent des 
informations relatives & des congestions dans le reseau. 

10. Dispositif selon l'une des revendications 7 h 9, carac- 
25 terise en ce que les secondes donnees comprennent des 

informations relatives a des depassements de deiais (TO) dans 
le reseau. 

11. Dispositif selon l'une des revendications precedentes, 
30 caracterise en ce que ladite matrice produit est representa 

tive d'un debit dans le reseau associe a la connexion & 
surveiller ou h simuler. 

12. Dispositif selon la revendication 11 , caracterise en ce 
35 que ladite matrice produit est representative d'un debit 

moyen dans le reseau. 
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13. Dispositif selon la revendication 11, caracterise en ce 
que ladite matrice produit est representative de fluctuations 
d'un debit instantane. 

5 14. Dispositif selon l'une des revendications pr6c<§dentes , 
caracterise en ce que le module de modeiisation est agence 
pour construire une plurality de matrices, en nombre corres 
pondant sensiblement au nombre de paquets dans le r§seau. 

10 15. Dispositif selon l'une des revendications pr6c6dentes, 
caracterise en ce que ledit module choisi comprend la 
consideration de la taille variable d'une fenetre (w n ) 
utilis^e pour contrdler le nombre de paquets dans le r6seau. 

15 16 - Dispositif selon la revendication 15, caracterise en ce 
que ledit modele choisi comprend la consideration d'un reseau 
comprenant des routeurs k discipline de type "premier arrive 
premier s ervi " . 

20 17- Dispositif selon l'une des revendications 15 et 16, 
caracterise en ce que ledit module choisi comprend la 
consideration d'un reseau comprenant des routeurs & disci- 
pline de type WFQ. 

25 18 - Dispositif selon l'une des revendications precedentes, 
caracterise en ce que le module choisi comprend la conside- 
ration d'un reseau controie par un protocole de type TCP. 

19. Dispositif selon l'une des revendications 15 & 18, 
30 caracterise en ce que ledit module choisi comprend la 

consideration d'un service (s) deterministe. 

20. Dispositif selon l'une des revendications 15 h. 18, 
caracterise en ce que ledit module choisi comprend la 

35 consideration d'un service (s) aieatoire. 
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21. Dispositif selon 1'une des revendications precedentes, 
caracterise en ce que la ma-trice courante au moins comporte 
des coefficients (A ik ) dynamiquement variables, 

5 22. Dispositif selon 1'une des revendications precedentes, 
caracterise en ce que la matrice produit obtenue est un 
vecteur repr^sente par une matrice & colonne unique (Z(n)). 

23. Procede d'aide a la surveillance d'un systeme complexe, 
notamment d'un reseau de communication , caracterise en ce 
qu'il comprend les etapes suivantes : 

a) obtenir des premieres donnees representatives de pararae- 
tres du reseau, 

b) construire une premiere matrice (Z(n)), selon un modele 
choisi, en fonction desdites premieres donnees , 

c) recevoir, a un instant choisi, au moins de secondes 
donnees representatives d'evenements dans le reseau, 

d) construire au moins une seconde matrice (A(n)) de struc- 
ture dynamiquement variable , selon le modele choisi, en 
fonction des secondes donnees, et 

e) effectuer sur lesdites matrices une operation formant 
produit selon l'algebre dite MAX-PLUS , la matrice produit 
obtenue (Z(n+1)) etant representative de l'etat du reseau 
audit instant choisi. 

24. Procede selon la revendication 23 , caracterise en ce 
qu'il comporte en outre l'etape suivante z 

f) repeter, a des instants choisis, les etapes c) , d) et e), 
tandis que la matrice produit obtenue devient la premiere 
matrice apres l'etape e) , 

ce qui permet de suivre une evolution temporelle de l'etat du 
reseau auxdits instants choisis. 

25. Procede de simulation d'un systeme complexe, notamment de 
35 mecanismes et de controles de flux et de congestions dans un 

reseau de communication, caracterise en ce qu'il comprend les 
etapes suivantes : 

a) obtenir des premieres donnees representatives de parame- 
tres propres au reseau, 
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b) construire une premiere matrice (Z(n)), selon un modele 
choisi, en fonction desdites premieres donnees, 

c) simuler des evenements dans le reseau et prevoir au moins 
de secondes donnees representatives desdits evenements, 

5 d) construire au moins une seconde matrice (A(n)) selon le 
modele choisi, en fonction desdites secondes donnees, et 

e) effectuer sur lesdites matrices une operation formant 
produit selon l'algebre dite MAX-PLUS, la matrice produit 
obtenue (Z(n+1)) <§tant representative d'un <§tat du reseau 

10 subissant lesdits Evenements. 

26. ProcEde selon la revendication 25, caracterisE en ce 
qu'il comporte en outre l'Etape suivante : 

f) repEter, pour des evenements successifs, les etapes c), d) 
15 et e ) , tandis que la matrice produit obtenue devient la 

premiere matrice aprEs l'Etape e), 

ce qui permet de prevoir une Evolution de l'Etat du rEseau en 
fonction desdits EvEnements. 
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